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Abstract 

We write down the gap equation for the fermion self-energy in a finite- 
temperature abelian gauge theory in three dimensions. The instantaneous ap- 
proximation is relaxed, momentum-dependent fermion and photon self-energies 
are considered, and the corresponding Schwinger-Dyson equation is solved nu- 
merically. The relation between the zero-momentum and zero-temperature 
fermion self-energy and the critical temperature T c , above which there is no 
dynamical mass generation, is then studied. We also investigate the effect 
which the number of fermion flavours Nf has on the results, and we give the 
phase diagram of the theory with respect to T and Nf. 
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1 INTRODUCTION 



In this paper we study dynamical fermion mass generation in a three-dimensional 
abelian gauge theory at finite temperature. The interest in phenomena associated 
with finite-temperature strongly coupled systems lies not only in their non-trivial 
field theoretical interpretation but also in the fact that in some cases specific compar- 
isons between experimental and theoretical results can be made. Dynamical fermion 
mass generation is such a phenomenon whose role in such a priori unrelated effects 
as chiral symmetry breaking in QCD or superconductivity in solid state systems 
gives it a particular importance. In the latter case, given the fact that certain cop- 
per oxides exhibit almost two-dimensional high-temperature superconductivity, the 
study of 2 + 1 field theories can be particularly instructive. 

Perturbation theory on its own can unfortunately do little in exploring the 
critical behavior of such systems. One has to use therefore non-perturbative tech- 
niques like effective potentials and Schwinger-Dyson (S-D) equations, which, al- 
though they cannot guarantee the precision of their quantitative results, provide 
interesting qualitative insights. 

We proceed therefore by writing down the S-D equation for the fermion self- 
energy in the real-time formalism and then study its behaviour with temperature. 
In particular, of considerable interest is the quantity r = 2S(0)/fcsT c , where S(0) 
is the fermion self-energy at zero temperature and momentum and T c is the critical 
temperature above which there is no dynamical mass generation. The identification 
of A = 2S(0) with the superconductor gap makes the model in principle experi- 
mentally testable Q. 

Similar investigations have been hitherto limited by the complexity of the 
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problem, in that the fermionic and photonic two-point functions which appear in 
the gap equation depend not on only one variable -as in the zero-temperature case- 
but on two independent ones. Responsible for that is the preferred reference frame 
(heat bath) associated with the temperature, which breaks Lorenz invariance. Thus 
the quantities of interest depend independently on energy and three-momentum. 

Popular approaches so far use drastic approximations such as taking the 
fermion self-energy to be momentum independent and truncating the energy de- 
pendence of the photon self-energies, either by considering them as functions only 
of the four-momentum squared Q , or by taking the energy to be zero (the "instan- 
taneous approximation") 

The study presented here is the first attempt to relax simultaneously all 
these approximations by numerically solving the problem and keeping the correct 
momentum dependence of the self-energies. This allows us to explore the critical be- 
haviour of the theory with respect to temperature and number of fermion flavours. 
However, due to the complexity of the equations, we neglect their imaginary parts, 
even though they could play an important role ||. Moreover, we consider only the 
one- fermion- loop contribution to the photon polarization diagrams, an approxima- 
tion justified for a large number of fermion flavours. We also choose to truncate the 
infinite S-D-equation hierarchy by replacing the full photon-fermion vertex by the 
bare one. Being aware of the severity of this approximation in relation to gauge- 
invariance and wave- function renormalization Q, we plan to make use of a more 
suitable truncation of the S-D hierarchy in a future publication. 
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2 THE GAP EQUATION 



The Lagrangian of the three-dimensional abelian gauge theory under study is a 
variation of the usual QED one with massless fermions: 

where, as usual, F^ v = d^A u — d u A^, a = 1, Nf with Nf the number of fermion 
flavours, the representation of the gamma matrices is four dimensional (a possible 
choice is 70 = diag(i<7 3 , — ia 3 ), 71 = diag(iai, —i<?\), 72 = diag(icr 2 , -^2), where crj 
are the usual Pauli matrices) and 

/ 1 \ 
1 
T3 " 0-10 
\ -1 j 

Dynamical mass generation in this case is parity conserving, which could be energet- 
ically preferred to a parity non-conserving one @. We will refer to the corresponding 
model as T3 — QED, as it has been used in Q in connection with high-temperature 
superconductors. 

The presence of the T3 matrix is responsible for the fact that this model 
does not have chiral symmetries that would be broken by a non-zero local order 
parameter like < ^l^r >, which is associated with mass generation. Therefore, 
dynamical mass generation in this three-dimensional theory can occur at finite non- 
zero temperatures without breaking any global symmetries and therefore problems 
with the Mermin- Wagner theorem H are avoided. On the other hand, this will not 
influence the form of the gap equation, which is the same as for usual QED. 

The corresponding S-D formalism for the two-point fermion function in the 
rest frame of the heat bath gives the following gap equation in the real-time for- 
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malism for the fermion self-energy S, for external momentum p^ = (po,p) and loop 
momentum k^ = (ko,k) in Euclidean space: 

S(po, \P\) = W J j^Dpfa, \k\)S p ( P0 - k , \p- k\) (2) 

where Dg and Sg are the photon and fermion propagators, with the subscript (5 
indicating their temperature dependence, and a = e 2 Nf is the dimensionful coupling 
of the superrenormalizable theory under study. It can be considered as an effective 
ultra-violet cut-off of the model ||. In order to simplify the notation we write 
the first argument of the various functions symbolically as po instead of the more 



accurate y — p^. Since we are working at the one-loop level, an approach justified 
for large Nf, we do not expect field-doubling problems associated with the real-time 



formalism to play any direct role [1C]. 



A bare photon-fermion vertex is used for simplicity, since a more involved 
vertex would lead us to a system of coupled integral equations. The propagators 
Dg and Sg in the Landau gauge at temperature T are given by 



P ±£, T k*+Ilp(koAk\) ePM-\ 

^ ' U \k 2 + E 2 (k ,\k\) e^ol + l J y 0,1 u v ; 

where k 2 = k^k^, (3 = 1/ksT, we sum over the longitudinal and transverse photon 
polarizations P = L,T, we have dropped the ^ term of the fermion propagator 
since it disappears after the momentum integration, and we have suppressed the 
imaginary i's since they appear on both sides of the gap equation. 



4 



The photon polarization functions Hl,t appearing above are given by Q 
^ /, ,r,x otk „ 

u L (ko,\k\) = ^- + n ifc 
n T (fc ,|&|) = — -n ifc + n 2fc ,with 

o 

2ak 2 f d\p\ ( ( B k + D k y/ 2 \ 
2a f d\p\ ( f k 2 (B k + D k ) \ 1/2 \ 

U2k = -J s^y-y-^Di — ) ) 

B k = k 2 -ip 2 

D k = {Bl + 16^ 2 ) 1/2 (4) 

where the subscript k is just a reminder that the quantities are momentum depen- 
dent. It is worth noting that these are calculated via a one-loop massless fermion 
diagram, using the fermion propagator of Eq.||. The quantities Ti\^K are the finite- 
temperature contributions to the photon polarization. They provide the thermal 
screening responsible for the softening of the infrared behaviour of the theory, since 
at small loop momenta they take values on the order of what is usually referred to 
as "plasmon mass" squared w| = ^ §. The one- loop approximation is justified 
for large Nf. However, the masslessness of the fermions in this calculation could 
introduce in principle consistency problems with Eqj2] which describes fermion mass 
generation. We will return to this issue later. 

We are thus confronted with a three-dimensional non-linear integral equation 
for a function of two variables. An analytical study of the full problem seems an 
impossible task, so one is led to a computer simulation . Before proceeding to 
a numerical solution however, we can try to somewhat simplify the equation by 
noting that the delta function appearing in the photon propagator gives a negligible 
contribution. In fact, this function has two roots, one at very large momenta -of 
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order a- where the fermion self-energy, a decreasing function with momentum, is 
vanishingly small, and one at very low momenta where the integrand is also small. 



Several authors make the approximation of dropping this delta-function 0],|12],[13], 



and we will also adopt it. The S-D equation therefore takes the following form 

a f dk \k\d\k\d9 E{p - k ,\p- k\) ^ 1 



NfJ (27T)3 ( p _ k) 2 + E2(po _ ko:W _ k]) p f^ T fc2 + n P (*b,|fc|) 

/• |fc|d|fc|rffl |p-fc|) y> l 

NfJ (2vr)2 2E(eP* + 1) * e Jl 1P ^ T (p - eEf + P + U P (p - eE, \k\) ^ 

where we sum over the photon polarizations P = L,T and over the two roots 
of the delta function by introducing e = 1,-1, and E is given by the relation 
E 2 = \p-k\ 2 + £ 2 (£, \p-k\). 

We note that the definition of E involves the fermion self-energy, which we 
are trying to solve for, calculated at po = E. Therefore, we make the approximation 



E w \J\p- k\ 2 + H 2 (0,0), which is valid in the limits [] \p - k\ 2 > or < E 2 (0,0). 
Another issue related to E is that, for small external momenta, we have to calculate 
the photon polarization functions at regions which, if continued to Minkowski space, 
would make these functions imaginary O]. This is due to the fact that, whereas 
we study fermion mass generation, the polarization functions were calculated for 
massless fermions. In order to use consistently massless fermions in connection 
with the photon polarization, we should really calculate the functions lip at energies 
Po±\p—k\ instead oipo±E. We checked that both prescriptions give similar results. 

On physical grounds, we expect that increasing temperature ought to have 

decorrelating effects, which, after a certain point, should make the order parameter 

x We assume of course that E(0, |fc|) ~ E(E(0, 0), |fe|), since in the infrared, non-perturbative, 
region our formalism does not have the required accuracy. We checked however that our final 
results do not depend on this detail. 



< V'lV'-R >> an d the fermion self-energy associated with it, vanish. In fact, we see 
that the effect of finite temperature, both directly via the negative second term of 
Eq.||, and indirectly via the photon polarization functions is to reduce the integrand 
for a given self-energy. It is therefore expected that, for large enough temperatures, 
the self-energy will not find enough support from the right-hand side and dynamical 
mass generation will be impossible, i.e. the only solution to the gap equation will 
be the trivial one. We proceed now to verify -and quantify- this expectation. 

3 NUMERICAL RESULTS 



In order to attack the problem numerically, we have to discretize our external and 
loop momentum space having in mind the characteristic energy scales of the model 
as well as the expected behaviour of the relevant quantities. The coupling a sets an 
effective UV-cut-off, as has been seen in previous zero-temperature studies @. At 
small loop momenta, two quantities compete for the role of a physical IR-cut-off: the 
fermion self-energy £ and /cbT, the energy associated with the finite temperature. 
The latter would become in effect only when both external and loop-momenta, 
as well as the fermion self-energy, are small. Problems related to the ambiguity 
of this cut-off due to the non-analyticities of the photon polarization functions are 



discussed in Ref.[2] as well as in Ref.[ll]. This does not influence our results since in 
our study the fermion self-energy finally prevails as an IR-cut-off, as we will see that 
convergence of the algorithm is lost before £ becomes less important than ksT. The 
gap equation is therefore both IR- and UV- finite, and we will use IR and UV cut-offs 
only for numerical reasons. Noting that £ is a function decreasing with momentum, 
and that the expected hierarchy of scales is large (three orders of magnitude or 
more), we discretize p%, \p\ 2 , k%, \k\ 2 according to log 10 (Aj R ) + ^ log 10 (h? uv /h? IR ) 



and the angle 9 according to 2Tri/n, where i = 1, ...,n. We do not take the points 
9 = 0, 7r as integration points, since the kernel has an integrable singularity there. 

We thus have a lattice with five dimensions, two coming from the external 
and three from the loop momenta. Furthermore, we are always careful that Ajr < 
IcbT, £(0, 0) (by E(0, 0) of course we mean here and in the following S(A//j, A/#)), 
and A.uv > a > an d check the dependence of our results on the particular values of 
the cut-offs and the size of the momentum lattice, which we took to be 8 5 , 12 5 and 
16 5 . Even though the quantities E(0, 0)/a and kBT c /a vary with different choices of 
cut-offs or relaxation-speed parameter, the ratio r = 2£(0, 0)//cbT c , remains pretty 
stable. 

The numerical relaxation method employed to solve the gap equation consists 
of inserting an initial "input" configuration for the fermion self-energy S to (|), tak- 
ing the "output" configuration as the new "input" , and then iterating the equation 
until it is satisfied to a good accuracy. In particular, we consider that the algorithm 
has converged when the mean difference as well as the standard deviation of the 
points of the "input" and "output" functions is less than 10%. This is about the 
best accuracy our algorithm can achieve at zero temperature and Nf = 1. It should 
be noted that in all cases studied, for sufficiently small values for the temperature 
and fermion flavors, the mean difference of the input and output configurations con- 
verges much faster to zero than the standard deviation between them. This allows 
us to distinguish the deviation caused by numerical error from the deviation caused 
by a possible overall fall of the solution towards zero. One can then be confident 
that the solution found is stable and not converging to the trivial one. In order to 
avoid convergence to unstable, oscillating solutions corresponding to small £(0,0) 
M, we take the original configuration to be flat and equal to a. Initial configura- 
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tions falling smoothly with momentum only decrease the convergence time and are 
therefore preferable, since they are closer to the final solution, and they give the 
same results. 

The fermion and photon self-energies are needed inside the integrand on 
points outside the original momentum lattice. Their corresponding values are 
found by linear interpolation from the values of the quantities on the lattice. For 
large internal and external momenta, some arguments can fall also outside the 
lattice space. Even though we use some extrapolation to recover the functions 
with such arguments, the values of the self-energies for momenta near the UV- 
cut-off, even though small and insignificant for the final results, should be taken 
with caution. Moreover, in order to obtain a smooth solution we take the val- 
ues of the output function at an external momentum lattice point to be the 
average of the integration result S out , i.e. = (S out (z + + S ou t(* — 1 , j) 

+S out (i, j + 1) + S out (i, j — 1)) /4. This is just one of several similar stabilizing pro- 
cedures one could use, and we have checked that the results do not depend on the 
particular choice of such a procedure. 

For T = and Nf = 2 we get the solution shown in Fig. 1. It exhibits 
roughly the expected behaviour T,(p) ~ — ^/^°^ +S ^ for X(0) < p < a, with 7 
and 5 in principle functions of Nf Q. We find that the general form of S(j>o> 
does not change with various T or Nf. 

When the temperature exceeds some critical value T c , the algorithm does 
not converge, i.e. the deviation of the "output" from the "input" function oscillates 
enormously from iteration to iteration, both in magnitude and in sign, instead of 
decreasing monotonically and reaching the required 10% value. For these high tem- 
peratures, even though the algorithm does not converge, the "output" configuration 
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Figure 1: The self-energy £(po> \P\) at T = and Nf = 2 for a 16 5 lattice, plotted 
as a function of momenta po and \p\ in logarithmic scale. All quantities are scaled 
by the coupling a. 

for the self-energy after each iteration is tending on the average rapidly to zero (the 
higher the temperature, the more rapid the fall). Stating this differently, only a 
trivial input configuration seems to be able to satisfy the self-consistency of the 
gap equation. On the other hand, for temperatures just below critical, £(0,0) is 
still at about the same the order of magnitude as its zero-temperature value. This 
does not necessarily mean that the self-energy drops "immediately to zero" when 
the critical temperature is reached. It most probably means that we are slightly 
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Fermion flavours =>■ 


N f = l 


N f = 2 


N f = 3 
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0.68 


0.11 
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16 5 


21 


4.8 


8.7 


2.4 


0.4 


12 


0.23 


0.035 


13.1 



Table 1: The quantities s = 10 3 x E(0,0)/a at T = 0, t c = 10 3 x k B T c /a and 
r = 2so/t c for different numbers of fermion flavours and lattice sizes n 5 and for 
A uv /a = 0.1. 
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23 


5.1 
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3.3 


0.54 


12.2 


0.5 


0.074 


13.5 



Table 2: The quantities s = 10 3 x S(0,0)/a at T = 0, t c = 10 3 x k B T c /a and 
r = 2so/t c for different values of flavours and ultra-violet cut-offs A;yy, for a 16 5 
lattice. 

underestimating T c , since at temperatures just below the real critical value it is very 
difficult to find an input configuration irregular enough to achieve convergence to 
the solution, which in the neighbourhood of T c is probably oscillating and on the 
average much smaller in magnitude. 

We list in Table 1 the values of the fermion self-energy £(0,0) at zero mo- 
mentum and temperature, of the critical temperature, both scaled by a, and of 
the ratio r, for different values of the number of fermions Nf and lattice sizes. We 
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choose the UV-cut-off such that Ajjy/a = 0.1, so that not only it is smaller than 
the physical cut-off a, but also truncates the integration before the uninteresting 
region where the self-energy is vanishingly small. Even though the X and T c values 
do not follow a specific pattern with respect to the lattice size n 5 , the r ratio seems 
to be converging for larger n. In particular, the difference between the r-values 
for the 12 5 and 16 5 lattice is consistent with the 10% accuracy provided by the 
numerical algorithm. However, given the large scale hierarchies of the problem, the 
n = 8 case for instance is just indicative. Moreover, since we justified the one-loop 
calculations with the 1/Nf expansion, the small values of Nf, and especially the 
ones reported for Nf = 1, should be taken with caution. The values for Nf = 2 
are the ones relevant to the high-T c superconductors. Furthermore, we see that for 
Nf = 3 the hierarchy between S(0,0) and a is very large. The relatively small 
lattice size forced upon us by computer-power limitations might therefore explain 
the stronger dependence of so and t c on the lattice size for three fermion flavors. 

In Table 2 we present our results for a 16 5 lattice for different ratios Ajjv l a - 
The ratio r follows a similar pattern for various UV-cut-offs, which is expected from 
a super-renormalizable theory. Lower values for Ajjy would truncate the integration 
too early, since the self-energy would be non-negligible there. When A;yy is taken 
to be larger than a, the algorithm does not converge, showing that it is not easy to 
perform the integrations on momenta above the physical cut-off a with such small 
lattices. The reason is that numerical errors dominate the calculation at this region, 
since the self-energy, a function decreasing with increasing momenta, is very small 
there (typically S(0, 0)/a is on the order of 10 -2 and smaller). 

We see that the ratio r increases with increasing Nf. This behavior, as well 
as the overall magnitude of r, is consistent with the results of Ref.[^] for the case 
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including retardation, when the photon polarization functions are approximated by 
II L ~ « ^ + 2t(jp or ^ + Wp, as in [Q. Since in our treatment the exact one- 
loop photon polarizations for massless fermions are used, our results indicate that 
this is the most sensible approximation. The other approximations in [Q] get smaller 
r ratios by increasing at the same time the value of the critical temperature by a 
factor as large as 2 or 3, which renders their validity questionable. An indication 
that the results in Tables 1 and 2 cannot be easily compared to the ones in |^] for 
various numbers of fermion flavors is the fact that in that paper for instance results 
for Nf = 4, 5 are reported, whereas in our case the behaviour of the theory with 
Nf renders mass generation with so many fermion flavors difficult or very small, 
which is in accordance with theoretical expectations. A more accurate study of the 
behaviour of the theory with Nf is enabled by the relaxing of several approximations 
used in Q, and is presented later in this paper. 

The values given for r could be overestimating this ratio, as a result of 
several effects. First, since £ is momentum dependent, one could take the relevant 
value for the ratio r to be S(p = S(0)) instead of S(0). This could decrease r by 
10 — 20%, which is justified not only by the accuracy of our algorithm but also by our 
limited understanding of infrared dynamics. Second, as was already noted, at the 
temperature at which fluctuations destroy the convergence of our algorithm, S(0, 0) 
has still the same order of magnitude as its value at zero-temperature. This not 
only shows that the fermion self-energy falls pretty abruptly when the temperature 
reaches its critical value, but that we might be slightly underestimating T c . The 
self-energy might be non-zero and falling with temperature for somewhat higher 
temperature values by taking highly irregular shapes that our algorithm is unable 
to find. In addition, it might not be fair to expect the algorithm at zero temperature 
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and at temperatures close to the critical value to reach the same accuracy of 10%, 
which after all determines which value of critical temperature we report, since larger 
temperatures make the numerical error coming from the second term of the right- 
hand side of Eq.|| more important. 
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Figure 2: The fermion self-energy at zero momentum and zero temperature, scaled 
by a and on a logarithmic scale, with respect to Nf for a 16 5 lattice. We fit our 
results with the curve £(0,0) /a = exp {-4.5/(N c /N f - l) 1/2 }/6, with N c = 4.35. 
Values of Nf larger than 3.4 are not considered, because then the self-energy falls 
below the IR-cut-off. 



We found that for large values of Nf the algorithm is also not converging. 
As in the critical temperature case, failure of convergence is accompanied with a 
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rapid decrease of the output configurations, which is a sign that the solution tends 
to the trivial one. This indicates that there might be a critical value iV c above which 
dynamical mass generation is impossible. In particular, for a lattice size of 16 5 we 
found that N c = 3.4 for k B T « E(0, 0) Air ps 10~ 5 x a. Larger values of Nf not 
only produce fermion self-energies tending below the IR-cut-off, but are also unable 
to reach the required accuracy of 10%. This value is consistent with theoretical 
expectations which give N c = 32 /it 2 at zero temperature ||. 

In Fig. 2 we plot the fermion self-energy at zero momentum and zero tem- 
perature as a function of Nf, and we fit it with the -phenomenological- curve 
E(0, 0)/a = exp {-A.5/(N C /N f - l) 1 / 2 }/6. From the fit a value N c = 4.35 seems to 
be close to our data. The non-analytic form of the fitted curve is needed not only 
to describe the exponential fall-off of the self-energy but also its probable vanishing 
at the critical value of Nf, and it exhibits the non-perturbative character of the 
dynamics. However, the approach of the self-energy to Air does not allow us to 
draw firm conclusions on the exact vanishing of £ at a critical number of fermion 
flavors. The deviation of this formula from theoretical expectations which give it as 
£(0,0) /a » exp{-2ir/(N c /N f - l) 1 / 2 } with N c = 32/ V 2 §§ is most probably due 
to the fact that the latter non-analytic formula is intended only for Nf very close 
to N c , and is not expected to fit data with smaller Nf. 

In Fig. 3 we proceed in a similar manner for a non-zero temperature case, 
namely kRT c /a = 10 -4 . We fit the zero- momentum fermion self-energy with the 
-phenomenological- curve E(0,0)/a = racp{-2.1/(JV c /JV/ - l) 1 / 2 }/^. From the fit 
a value N c = 3.2 is favored. However, we loose convergence of our algorithm at 
Nf ~ 2.6, which should be much closer to the real critical value -/V c , according to 
the criteria used so far. This deviation indicates, as in the zero-temperature case, 
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that the non-analytic functional form of the fit should not be used for values of 
Nf far away from N c . In this finite-temperature case, convergence of the numerical 
algorithm is lost before the self-energy falls below the IR-cut-off. Figures 2 and 3 are 
similar to the corresponding ones in [||, which also indicate the existence of a critical 
behaviour with Nf, even though our value for N c at zero temperature is closer to 
the theoretical expectations. The decrease of N c with increasing temperature is 
also expected, and is shown in more detail in the following figure, i.e. the phase 
diagram. 

In Fig. 4 we plot the phase diagram of the theory with respect to Nf and 
JzbT, where by the term "phase" we mean a situation where there is or there is no 
dynamical mass generation. The non-perturbative nature of our formalism allows 
us to consider also non-integer values of Nf . We were able to fit the critical line with 
the -phenomenological- curve fcs^c/a = exp {— 2.5-/Vy}/17. A similar exponential 
fall-off has already been seen in [||. This functional form should not be used for 
values of Nf larger than about 3, because then T c falls below the IR-cut-off. It is 
quite possible that when Nf ~ 3.4 the behavior of ksT with Nf becomes here also 
non-analytic. Nevertheless, we are not able to explore this region, since for values 
of Nf larger than about 3, kT c falls below the IR-cut-off. 

4 CONCLUSIONS 



We report results indicating that for temperatures below a finite critical value, 
dynamical mass generation in finite temperature T3 — QED in three dimensions is 
possible. We find a large ratio r that not only exceeds the usual BCS prediction, 
but also the value measured in high-temperature superconductors, which is close 



to 8 [14], and indicates anyway that these are very strongly coupled systems. This 
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ratio, on the order of 10, is generally consistent with results of previous studies [0], 
even though its precise behaviour with Nf shows some differences. We also find 
that for a number of fermion flavours larger than a critical value, of order 3, there 
cannot be mass generation, or, even if there is any, it falls below any IR-cut-off 
that we may set so it is impossible to study with the present algorithm. We are 
then able to draw the phase diagram of the theory, with respect to temperature and 
fermion flavours, which separates the regions where there is and where there is not 
dynamical mass generation. 

We were able to go beyond the instantaneous approximation and consider 
momentum dependent fermion and photon self-energies by working with quantities 
depending on two independent variables simultaneously. Previous results on this 
context, relying on severe truncations and approximations, are thus put now on a 
firmer basis. The increased lattice dimensionality that followed however allowed us 
to use only a limited number of lattice cites in each dimension, even though larger 
lattices would be better suited for a problem with such large scale hierarchies as the 
present one, and they would possibly allow for a better accuracy than 10% in our 
solutions. 

It is also true that the specific truncation of the S—D equation that we chose, 
as well as the dropping of the imaginary parts of all the self-energies could influence 
these results non-trivially. Moreover, the number of fermion flavours might not be 
large enough to justify the 1/Nf expansion and the one-loop diagrams we consider. 
We will try to return to these issues in a future publication. 
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Figure 3: The self-energy at zero momentum and at ksT/a = 10 4 , scaled by a 
and on a logarithmic scale, with respect to Nf for a 16 5 lattice. We fit our results 
with the curve S(0,0)/a = exp{-2.1/(JV c /JV/ - l) 1/2 }/15, with N c = 3.2. The 
algorithm does not converge for values of Nf larger than 2.6. 
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Figure 4: The phase diagram of the theory with respect to ksT on a logarithmic 
scale and Nf for a 16 5 lattice. We fit the critical line with the curve ksT/a = 
exp {— 2.5iV/}/17. This functional form could lose its validity for Nf larger than 
about 3, since then T c falls below the IR-cut-off. 
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